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Anatomy of A Dark Burst - The Afterglow of GRB 060108 
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£3 . ABSTRACT 

We present a multiwavelength study of GRB 060108 - the 100th Gamma Ray Burst 
discovered by Swift. The X-ray flux and light curve (3-segments plus a flare) detected 
with the XRT are typical of Swift long bursts. 

We report the discovery of a faint optical afterglow detected in deep BVRi' band 
imaging obtained with the Faulkes Telescope North (FTN) beginning 2.75 minutes 
' after the burst. The afterglow is below the detection limit of the UVOT within 100s of 

the burst, while is evident in A'-band images taken with the United Kingdom Infrared 
Telescope (UKIRT) 45 minutes after the burst. The optical light curve is sparsely 
sampled. Observations taken in the R and i' bands can either be fit with a single 
power law decay in flux, F(t)oc t~ a where a = 0.43 ± 0.08, or a 2-segment light 
curve with an initial steep decay ai<0.88±0.2, flattening to a slope a2~0.31±0.12. A 
marginal evidence for rebrightening is seen in the i' band. 

Deep i?-band imaging obtained ~ 12 days post burst with the VLT reveals a 
faint, extended object (R ~ 23.5 mag) at the location of the afterglow. Although 
the brightness is compatible with the extrapolation of the slow decay with index a% 
significant flux is likely due to a host galaxy. This implies that the optical light curve 
had a break before 12 days, akin to what observed in the X-rays. 

We derive the maximum photometric redshift z < 3.2 for GRB 060108. We find 
that the Spectral Energy Distribution at 1000 s after the burst, from the optical to 
the X-ray range, is best fit by a simple power law, F„ oc v~P ', with fiox — 0.54 
and a small amount of extinction. The optical to X-ray spectral index (Pox) confirm 
GRB 060108 to be one of the optically darkest bursts detected. Our observations rule 
out a high redshift as the reason for the optical faintness of GRB 060108. We conclude 
that a more likely explanation is a combination of an intrinsic optical faintness of the 
burst, an hard optical to X-ray spectrum and a moderate amount of extinction in the 
host galaxy. 
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1 INTRODUCTION 

Gamma ray bursts (GRBs) are brief, intense and totally 
unpredictable flashes of gamma rays on the sky that are 
thought to be produced during the core collapse of mas- 
sive stars (long-duration bursts) or the merger of two com- 
pact objects such as two neutron stars or a neutron star 
and stellar-mass black hole. Until the recent launch of the 
Swift satellite in November 2004, it was notoriously dif- 
ficult to observe GRBs at other wavelengths within sec- 
onds or minutes after the burst. Nevertheless, the success- 
ful identification with BeppoSAX of br ight, long-lived X - 
ray afterglow emission for long bursts ijCosta et alJ I199t1) 
and th at of the correspon ding optical and infrared counter- 
parts llvan Paradiislll997ft . established GRBs as cosmolog- 
ical, and therefore the most instantaneously luminous ob- 
jects in the Universe. Similar breakthroughs for short bursts 
have recently occurred, showing them also to be extragalac- 
tic, but less luminous and less distant than long bursts 
iGehrels^^l . 2005 : IVillasenor et alJl2005t iFox et af]l2005l 
ICovino et alJl2006t iBarthelmv et al]l2005ft . 

With the availability of Swift's promptly-disseminated 
arcsec localizations and the on-board rapid-s lew X-ray and 
ultra violet / optical telescopes (XRT, UVOT; iGehrels et alJ 
2004), multi- wavelength monitoring of GRBs from the ear- 
liest possible times is now being performed for a significant 
number of bursts. Additionally, large aperture ground-based 
robotic telescopes such as the 2-m Liverpool (Steele 2004) 
and Faulkes telescopes respond rapidly to GRB alerts and 
begin automatically imaging the target field within minutes 
of receipt of an alert, providing early deep upper limits or 
multi-colour follow-ups of optical counterparts as faint as 
R~ 18 -22 mag (e.g. Guidorzi et al. 2005a, 2006b; Monfar- 
dini et al. 2006b). 

Despite increasingly rapid responses that provide sen- 
sitive limits within minutes of the burst, the absence of 
long-wavelength emission afterglows for a significant number 
of GRBs (so-called "dark bursts") remains a key unsolved 
problem. In the pre-Swift era, as many a s 50% of BeppoSax 
bursts were lacking an optical detection iDe Pasauale et all 
l2003l:lLazzati. Covino fc Ghiserlinl2002ft . The discovery rate 
of opt ical afterglows wa s even higher for HETE2 than for 
SAX jLamb et al.ll2004ft . and it was expected to increase 
significantly in the Swift era of rapid followup. Instead, 
a substantial fractio n of Swift bursts r emain undetected 
in the optical band jRoming et alll2006ft . Possible scenar- 
ios to explain the "observed" optical darkness of these 
bursts, apart from fast-fading t ransients lacking su fficiently 
deep, early-time observations llGroot e t al.||l9 98ft. include 
intrinsically-faint optical afterglows (|FaIb^-t_alJ [2001]; 
Lazzati. Covino fc Ghisellinill2002t iDe Pasauale et alJl2003: 
Berger et al 120051) . lighly obscured afterglows whose optical 



light is absorbed by the circumburst o r interstellar mate- 
rial llLazzati. Covino fc Ghiselliml 120021: IDe Pasauale et alJ 



20061). high redshifts (lFruchterll999l 


Lamb & Reichartl2000l 


Bromm & Loebl 120021: iBerger et al. 


120021 iTagliaferri et alJ 


2005), and radiative suppression in a sub-class of bursts with 



unusually high 7-ray effi ciency producing intrinsically low X- 
ray a nd optical fluxes jRoming et aljl2006t iPedersen et alJ 

Here we present a multi-wavelength X-ray, optical and 
infrared study of the optically-faint GRB 060108. The Burst 



Alert Telescope f BAT. IBarthelmv et al.ll2005ft was triggered 
by this GRB - Swift's 100th burst - at 14:39:11.76 UT on 
January 8 th 2006. The 7-ray light cur ve has a single pea ked 
structure with a FRED time profile dOates et al] 12006ft . a 
duration T 90 = 14.4 ± 1.6 s (in the 15 - 350 keV band), and 
a 15-150keV fluence S 7 = (3.7 ± 0.4) x 10" 7 ergcm -2 . 
Here and in the following errors are at 90% confidence 
level, unless specified otherwise. The Swift X-ray and ul- 
traviolet/optical telescopes (XRT/UVOT) began observing 
at 91 s and 76 s after the BAT trigger respectively, followed 
shortly after by the Faulkes Telescope North (FTN) at 2.75 
minutes post trigger. Infrared observations with the United 
Kingdom Infrared Telescope (UKIRT) were acquired dur- 
ing the optical imaging period beginning at 45 min post- 
burst. The Very Large Telescope (VLT) was used to ob- 
tain further near infrared imaging at 16.1 hours and 2.7 
days. Following an initial estimat e of the position of the 
X-ray a fterglow iPage e t al||2006ft. a revised location was 
derivedj|Bjrtl£r^^B^n^^Sj); the location of the op- 
tical |M onfta^mT^t^lT T2006air and infra red counterparts 
JP'Avanzo et alJl2006t iLevan et al"ll2006aft was found to be 
consistent with this revised XRT position. Deep _R-band 
imaging and spectroscopy were subsequently performed with 
the VLT at 12.7 days and ~21 days respectively. 



2 OBSERVATIONS AND ANALYSIS 

2.1 Swift Gamma- Ray Observations 

At 14:39:11.76 UT, on 8 th January 2006, the Swift BAT 
triggered and loca ted GRB 060108 (BAT trigger 176453) 
JOates et al]l2006ft . Unless otherwise specified, times in this 
section are referenced to the BAT trigger time, hereafter 
designated To- The burst was detected in the fully coded 
part of the BAT field of view, meaning that all BAT detec- 
tors were illuminated by the source. The spacecraft began 
to slew to the source location 11.2 s after the trigger and 
was settled at the source location at To+73.20 s. 

The BAT data for GRB 060108 between T -300s and 
To+300 s were collected in event mode with 100 /is time res- 
olution and 6keV energy resolution. The data were pro- 
cessed using standard Swift-BAT analysis tools and the 
spectra were fit using XSPEC 11.3. Before processing, the 
BAT event data were corrected for mask tagged weighting 
using the standard BAT tool, batmaskwtevt. This procedure 
scales each count in each detector by an amount propor- 
tional to the portion of the detector exposed to the source 
through the empty spaces in the BAT coded mask. Mask 
weighting effectively removes the background and all BAT 
light curves in this paper have been background subtracted 
by this method. 

2.2 Swift X-ray Observations 

The X-ray Telescope (XRT, iBurrows et"afl l2005aft po- 
sition refined analysis located thi s burst at q(J2 000) = 
09 h 48 m 01!6, 5(J2000) =31°55'04'.'6 iPage et al]l2006ft . XRT 
began observations approximately 90s after the burst trig- 
ger. The analysis of the XRT data was performed using the 
XRT pipeline software version 2.3. The initial exposure was 
taken in Image Mode and no source was detected. The XRT 
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then observed in Windowed Timing (WT) mode for ~ 10 s 
before switching to Photon Counting (PC) mode for a fur- 
ther 10 ks; a fading uncatalogued source was detected. Swift 
observed the field until 10 6 s after the burst trigger. 

The source and background counts used to create the 
light curve and the spectra were obtained from the cleaned 
event files using extraction regions. The analysis of the WT 
mode data provided only limited results as the exposure was 
of such a short duration. 

The first 1000 s of PC data were found to be piled up. 
The size of the region that was affected by the pile-up was 
determined by comparing the measured point spread func- 
tion (PSF) with a model XRT PSF. The affected radius was 
found to be 1.6 pixels (4") in radius, therefore an annular re- 
gion with inner radius 4" and outer radius of 30 pixels (70.8") 
has been used for the extraction. The rest of the PC data 
were extracted using a circular region of radius 30 pixels 
(70.8") and a background region of radius 80 pixels (189"). 
These regions were used to extract data from all orbits of 
PC data. 

When extracting the PC mode spectrum, only events 
with grade 0-12 were used, and for the WT mode, events 
were used of grade 0-2. The energy response of the detector 
was taken into account when fitting the spectrum. The re- 
sponse matrices (RMs) were taken from the Swift calibration 
database, CALDB 20060426. These were the most recent 
calibration products available and were released in April 
2006. The effective area files were created with xrtmkarf 
and required an image file created using expomap. This code 
maps the detector plane so any bad pixels or hot columns are 
accounted for in the effective area file. A correction factor 
calculated with xrtmkarf was applied to the piled up section 
of the XRT light curve. 

During the observations, the burst's location near to the 
Sun caused the XRT to warm up, leading to mode switching 
of the instrument. This only affected data in the second orbit 
and onwards. Thus the PC light curves, where appropriate, 
were also corrected for low fractional exposure. 



Table 1. Swift UVOT observations of GRB 060108, corrected for 
Galactic extinction. E(B - V) = 0.018. 
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Time Range 

to 
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to 
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(mag) 
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75-86 
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17.9 
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90-290 


200 


19.6 
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76-1096 
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16835-17735 


900 


20.8 


B 
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20.1 


B 


10983-40489 


1588 
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U 
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20.2 
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Figure 1. Faulkes Telescope North i'-band image of the optical 
counterpart to GRB 060108 (indica ted with an arrow). Th e orig- 
inal JPage et al.ll2006ft and revised feutler fc Bloomll2006T) XRT 
error circles are also shown. 



2.3 Optical and Near Infrared Imaging 

2.3.1 Swift Ultraviolet Optical Telescope (UVOT) 

The UVOT began observations 75 s after the BAT trigger. 
The first ~ lis exposure, taken while the spacecraft was set- 
tling on the target, was taken through the V filter in photon- 
counting ("Event") mode. Once the pointing had stabilized, 
a 200 s V-band finding-chart exposure in Image/Event mode 
was taken. After this, the UVOT cycled through the colour 
filters. The later exposures were in image mode only. Counts 
were extracted from the images using a 3" radius aperture at 
the position of the optical afterglow given in this paper. The 
count rates were aperture-corrected to ensure compat ibility 
with the photometry calibration bv lPoole et alJ <|2005T) . We 
then derived the 3<r upper limits in count units, which were 
converted to magnitudes using the zero-points available in 
the Swift CALDB, and corrected for Galactic extinction us- 
ing E(B - V) = 0.018. Results are given in Table 



2.3.2 Faulkes Telescope North (FTN) 

The Faulkes Telescope North, located on Maui, Hawaii, re- 
sponded robotically to the Swift alert and began observ- 
ing 2.75 min after the GRB trigger. The initial detection 
mode images (three 10 s i?-band exposures) were processed 
automatically by the Liverpool GRB pipeline, LT-TRAP 
jGuidorzi et al.ll2006aD . No new source has been identi- 
fied in the XRT error box dPaee et alJl2006l) brighter than 
R ~ 19.5 mag (At = 3.1 min) and R ~ 20.2 mag (At = 11.7 
min) in co mparison with USNO-B 1.0, 2MASS and GSC 2.3 
catalogues IIGuidorzi et alj|2006bf) . The telescope continued 
in multicolour imaging mode, obtaining the exposures in 
BVRi' bands. 

The robotic followup ended after one hour. Subse- 
quent post-processing of the i! band images revealed an 
uncatalogued source ~3" NE of the edge of the XRT er- 
ror circle (see Figure 0. The i'-band brightness of this 
source in the FTN image relative to the neighbouring 
object was greater than in the pre-burst SDSS z'-band 
image, suggesting it as the candidate afterglow. Subse- 
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Table 2. Optical detections and lower limits. Errors are at 90% confidence level. 



Mean Date 


Filter 


Telescope 


Exposure 


Seeing 


Magnitude 


Time Range 


At 


CUT) 






Time (s) 


(") 




(s) 


(s) 


Tan 8 fil ^77fi 


iDCbbCll-iD 


FTN 


inn 


1 12 


^> -Ly . yu 


269- 


-279 


97 A 


Jan 8.624165 


Bessell-B 


FTN 


390.0 


1.54 


22.31±0.26 


450- 


-1552 


1172 


Jan 8.643794 


Bessell-B 


FTN 


540.0 


1.74 


22.72±0.35 


2198- 


-3310 


2868 


Jan 8.614304 


Bessell-V 


FTN 


10.0 


1.29 


21.7±1.2+ 


315- 


-325 


320 


Jan 8.612498 


Bessell-R 


FTN 


10.0 


1.09 


>20.20 


159- 


-169 


164 


Jan 8.612996 


Bessell-R 


FTN 


10.0 


1.09 


>20.20 


202- 


-212 


207 


Jan 8.616757 


Besscll-R 


FTN 


30.0 


1.18 


21.23±0.44 


517- 


-547 


532 


Jan 8.620774 


Bessell-R 


FTN 


240.0 


1.47 


21.84±0.28 


164- 


-1186 


879 


Jan 8.633794 


Bessell-R 


FTN 


300.0 


1.55 


21.88±0.31 


1767- 


-2359 


2004 


Jan 8.649744 


Bessell-R 


FTN 


420.0 


1.74 


21.98±0.31 


2944- 


-3711 


3382 


Jan 8.614940 


SDSS-i' 


FTN 


10.0 


1.12 


20.9±0.7 t 


370- 


-380 


375 


Jan 8.619292 


SDSS-i' 


FTN 


100.0 


1.34 


21.69±0.38 


375- 


-891 


751 


Jan 8.630704 


SDSS-i' 


FTN 


300.0 


1.51 


21.45±0.16 


1352- 


-1993 


1737 


Jan 8.644338 


SDSS-i' 


FTN 


300.0 


1.56 


21.88±0.33 


2464- 


-3265 


2915 



t Marginal detection with low statistical significance. Errors are derived by adding statistical and 
systematic errors in quadrature. Magnitudes are uncorrected for Galactic extinction. 3<r upper 
limits are quoted for non-detections. Column 7 lists the Time range covered by the exposures and 
column 8 lists the corresponding exposure-weighted mean time (At) since the burst for co-added 
images calculated as At = tjA*i)/(^ At t ). 



quent analysis of the XRT data bv lButler fc Blooml tOOdh 
confirmed an XRT offset with a revised position for the 
afterglow a(J2000) = 09 h 48 m 0rl92, <5(J2000) = 31°55'07'.'8 
with 0'.'9 uncertainty (90% containment). As can be seen 
in Figure Q the location of the optical afterglow at 
a(J2000) =09 h 48 m 01 , !98, (5(J20 00) = 31°55'08'.'6 with ~ 0'.'5 
uncertainty (90% containment. IMonfardini et aljl2006bl) . is 
consistent with the revised XRT position. A similar lo- 
cation is derived for the NIR counterpart iden tified by 
iD'Avanzo et alJ j200eT) and iLevan et all J2006al) . Careful 
analysis of the optical data subsequently confirmed that the 
source was fading and further imaging with the FTN and 
other telescopes was manually triggered to follow the later- 
time evolution. Unfortunately moon and weather constraints 
prevented deep detections or limits beyond R ~ 23 mag. 

Table [5] summarizes the results of the optical and NIR 
imaging; as no FTN observations of standard stars were 
available, magnitudes in B, V, R, and i' were calibrated 
as follows. Observations in i' band were calibrated directly 
using the SDSS DR4 photometry 1 , while B, V, and R 
band data were cali brated assuming th e pre-burst SDSS 
g'r'i' photometry of ICool et alJ <l2006e 13) combined with 



Smith et al] Eo6 



observed filter transformations given by 
in their Table 7. Finally, the data were corrected for the 
Galactic extinction, which is low towards GRB 060108: 



A s = 0.09 mag, A v = 0.07 ma 


5 


Ar = 0.05 map 


A, = 0.04 mas dSchleeel et alJ Il99^ 




Cardelli et allll989l) 



Magni tudes were co nver ted into flux densities F (/iJy) fol- 
lowing [gessej| (|T979) and lFukueita et ail (I1996T) for compar- 
ison with X-ray flux densities. 



2.3.3 Very Large Telescope (VLT) 

J- and ifs-band observations of the field of GRB 060108 
were performed by using the ISAAC camera at two differ- 
ent epochs, and in the R band by using the FORS 1 in- 
strument at one epoch (see Table [3J. Both instruments are 
mounted on the ESO-VLT located at Cerro Paranal (Chile). 
All nights were photometric, with a seeing in the range 0.8"- 
1.4". Image reduction was performed following the standard 
procedures. Astrometric solutions were computed using the 
USNO catalog 2 . Aperture photometry was performed us- 
ing SExtractor feertin fc Arnoutslll996ri . In the optical, the 
flux calibration was performed observing L andolt standard 
stars . In the NIR, the 2MASS catalogue iSkrutskie et alJ 
2006) was adopted as a reference, even if only a handful 
of usable stars were present in the ISAAC field. Late-time 
low-resolution spectra were acquired with the VLT in an 
attempt to measure the host galaxy redshift. However, the 
object was quite faint and spectra have a very low signal-to- 
noise ratio, with no clear emission or absorption lines. We 
will not discuss further these spectra. 

In Table |H] we present the observing log and list the 
results of our spectrophotometric analysis for each epoch. 



3 RESULTS 

3.1 GRB 060108 and Its Afterglow 

The 



afterglow o f GRB 60108 is detected at X-ray 
JPaee et al) 120061) . optical (IMonfardini etHI l2006al) and 

near infrared wavelengths jLevaiiTt^n]2006a ) . The X-ray 
light curve shows the canonical 3-segment decay similar to 
many gamma ray bursts detected with Swift, with a small 
flare 300 s after the burst. Optically, the burst is not detected 



www .sdss.org/dr4/ 



2 http:/ /www. nofs.navy.mil/data/fchpix/ 
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Table 3. GRB 060108: VLT observation log. Errors are at 90% confidence level, while 3<r upper 
limits have been used for the non-detection. 



Mean UT 


Time since 


Exp 


Seeing 


Instrument 


Filter/ 


Magnitude 




burst (d) 


(s) 


(") 




grism 




Jan 9.307 


00.696 


10x90 


1.1 


ISAAC 


J 
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Jan 11.283 


02.672 
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0.8 


ISAAC 


Ks 


20.49± 0.15 


Jan 21.286 
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1.4 


FORS1 


R 


23.46±0.15 


Jan 29.237 


20.626 


3x2160 


0.8 


FORS2 


300V 




Jan 30.176 


21.565 


3x2160 


0.8 


FORS2 


300V 






5 10 
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Figure 2. BAT light curve of GRB 060108 in different energy 
bands: 15-25, 25-50, 50-100, 100-350keV, and total range 15- 
350 keV. 



at ultraviolet and optical wavelengths with the UVOT in the 
first two minutes, but it is detected in deeper B, V, R and 
i! band FTN images as early as 5.3 min after the burst. The 
optical afterglow of GRB 060108 is one of the faintest yet 
discovered. 



3.2 Swift Gamma-Ray Light curve and Spectra 

The BAT light curve (see Fig. |2J shows a simple triangu- 
lar profile in all energy bands, with a rise from background 
starting at Tb-4.0s, a peak at To, and a return to back- 
ground levels by Tb + 13.2s. A small secondary peak is seen 
at around To+10 s, but then no detectable emission has been 
registered after To+13.2s nor at the times of the flare seen 
by the XRT. The peak count rate measured by BAT is 2000 
counts s" 1 at To, and the duration is T90 = 14.4 ± 1.6 s 
(15-350keV band; estimated error including systematics) . 
The Is peak photon flux measured from Tb _ 0.48 s is 0.7 ± 0.1 
phcm" 2 s _1 (in the 15- 150 keV band). 

For the spectral analysis, we have analyzed the data 



between 15 and 150 keV, where BAT has a better response. 
We find that the emission from GRB 060108 shows the 
hard to soft spectral evolution commonly observed in many 
GRBs and that, for all time periods, the best fit to the BAT 
data is a simple power law. The spectrum of the first peak 
(Tb-4.02s to To+9s) is best fitted with a photon index of 
1.9 ±0.17 (x 2 / d -°- f - =54/90), while that of the smaller sec- 
ond peak (To+9s to Tb+13.2s) has a photon index 2.55 ± 
0.61 (x7 d -°- f - = 64/90). With the given spectral parame- 
ters, this burst appears to share some similarity with the 
X-ray flashes (XRFs). Although Swift is not well suited to 
identify XRFs, since it lacks coverage at low energy (i.e. in 
the the "classical" HETE band 2-30 keV), we notice that, 
if the spectrum is extrapolated down to this band, then the 
soft photon index would imply an XRF or at least an X-ray 
rich event. 



3.3 Swift X-ray Light curve and Spectra 

By following Ip'Brien et all fcOOd; . the BAT light curve can 
be extrapolated at lower energies and joined to the XRT 
light curve in a single plot. Although some consideration is 
needed (as the energy bands of the two instruments do not 
overlap), this is a powerful method which allows us to visual- 
ize the decay rate of the event through a longer time interval, 
commencing at the burst trigger. The combined BAT and 
XRT light curve resulting from this process is shown in Fig- 
ure |5] For the extrapolation, we used a BAT photon index 
of 1.9 (computed in the 15-150 keV band). 

The light curve has a complex shape, showing a steep 
decay, followed by a flatter phase which eventually steep- 
ens again. A visual inspection shows a possible flare, su- 
perimposed between the two segments. We modeled the X- 
ray light curve with both a double broken power-law and 
a double broken power-law with the inclusion of a Gaus- 
sian, to represent the flare. In the first case we find a 
X 2 /d.o.f. = 26/23, while by using the second model gives 
X 2 /d.o.f = 11/20. The probability that this improvement is 
due to chance, calculated using the F-test, is 0.05%. In or- 
der to check for overfitting effects, we repeated the fitting by 
adopting different choices for the temporal binning. We find 
that, in all cases, the light curve is better fitted by including 
a Gaussian, although in some cases both fits (with and with- 
out Gaussian) have reduced x 2 values larger than unit, and 
in turn the decrease in \ 2 ls l ess pronunced. The probability 
that the improvement is due to chance, calculated using the 
F-test, varies between 0.019 and 0.035 in the different fits, 
which means that the fit including the Gaussian can be used 
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Figure 3. BAT and XRT combined unabsorbed flux light curve 
of GRB 060108. 



Table 4. Best fit parameters for a fit of the XRT light curve 
of GRB 060108 with a double broken power law model plus a 
Gaussian. The epoch of the flare, tfl arc , and its width oq have 
been reported. From the Gaussian component, the Full Width 
Half Maximum is related to cfq by FWHM= 2.35cr(j. Errors are 
at 90% confidence level. 



Parameter 





2.78 ± 0.43 


CY'2 


0.30 ±0.07 


OC3 


1.05 ±0.10 


^break, 1 


250 ± 30 s 


^break,2 


(11.5 ±0.4) x 10 3 s 


^flarc 


300 ± 10 s 




32.3 ± 10.0 s 



with a confidence of ^ 2a. The best fit parameters appear 
to be always consistent with each others, and are given in 
Table [I] here and in the following we assume the convention 

F oc t- a v~ 13 . 

Our analysis shows that there are three distinct seg- 
ments in the X-ray light curve. The X-ray afterglow starts 
with a phase of steep decay (first segment, slope ~ 2.8), fol- 
lowed by a shallow section (second segment, slope ~ 0.3), 
and finally enters a "normal afterglow" decay phase, with 
a slope ~ 1 (third segment). The two break times which 
separate these three segments are tbroak.i ~ 250 s and 
tbreak,2 ~ 11-5 ks. This behavi our is typical of X- ray af- 
terglow lightcurves as shown bv lNousek et alJ <|2006T) . 

In order to perform a spectral analysis, we fitted the 
spectra taken during different time intervals by using an 
absorbed power-law model. The Galactic value of the col- 
umn density has been fix ed at N§ al = 1.7 x 10 20 cm~ 2 
.Dickey fc Lockmanlll990l this value is in agreement with 
that obtained fr om the more recent LAB survey, which is 
1.8 x 10 20 cm~ 2 . lKalberla et"alll2005l 3 ), and an additional 
absorption component, Nh, has been allowed to mimic the 
intrinsic excess in column density. Results are summarized 
in Tabled 

The spectral indices before and after the first break at 
250 ± 30 s have been first computed by including the flare 

3 http: / /www. astro. uni-bonn.de/~webrai/english/tools_labsurvey. 



at 300 ± 10s (Pi in the table). In the attempt to assess the 
flare contamination on these values, we divided the interval 
in three parts by isolating the time interval of the flare. Since 
these new spectra contained less than 100 background- 
subtracted counts, Cash statistics were used for the fitting 
(a minimum of 15-20 counts per bin are required for Gaus- 
sian x 2 statistics). A comparison of the resulting indices, P[, 
with the Pi previously computed, shows that the flare does 
not alter significantly the spectral indices registered during 
the first two decay phases. 

Because of the low statistics, when considering the spec- 
tra of individual segments separately we are unable to con- 
strain the value of the excess of column density over N§ al . 
We have only been able to derive a 90% upper limit, and 
only by using the spectra of the last two time intervals. In 
order to have an estimate for the photoelectic absorption, 
we re-extracted the X-ray spectrum using all data after the 
flare through to the end of the observations (time interval 
570 s< t < 4.63 x 10 s s). Since this spectrum has a bet- 
ter signal to noise, the excess Nh is actually significant at 
99% confidence level; the best-fit spectral parameters are 
Pi = l.06±° %, N H = 1-02± ;66 x 1()21 cm ~ 2 ( see also 

3.4 Optical and Infrared Properties 

The optical and infrared properties of the GRB 060108 af- 
terglow are also particularly interesting. Figure^shows the 
identification of the optical afterglow with the X-ray coun- 
terpart in an FTN i'-band image, while Figure 2] shows the 
optical and infrared light curves including detections and 
upper limits from optical to near infrared wavelengths, cor- 
rected for Galactic extinction. As shown in Table [5] the af- 
terglow is visible in all BVRi' bands from as early as At ~ 5 
minutes, but is not detected by UVOT within the first few 
minutes of the trigger or during the short exposure FTN R- 
band and B-band images taken at 164 s and 274 s after the 
trigger, respectively. Although these FTN exposures have 
only a ~ 10 s duration, the deep limits rule out a significant 
optical flare coincident in time with that possibly seen in 
the X-ray light curve. Moreover, the marginal FTN U-band 
detection (Table [5] and Figure [IJ reveals that the afterglow 
was already fainter than the UVOT limiting magnitudes. 

The optical light curves to 4000 s after the burst can be 
fit with a single power law decay in flux density F(t) oc t~ a 
where a = 0.43 ± 0.08 (Figure [5J, a slope that is also 
consistent with the gradient between the two A"-band de- 
tections. Alternatively, guided by the structure in the X-ray 
light curve and concentrating on the R and i' band light 
curves, a more complicated interpretation is possible. By a 
direct comparison between the X-ray and optical flux den- 
sity evolution (see Figure we find that, although the op- 
tical light curves are sparsely sampled (due to the need to 
coadd frames), there is evidence for a 2-segment light curve. 
A possible rebrightening is seen in i' band, although the 
significance is low and a monotonic decay is not ruled out. 
Using the early time i?-band upper limit as a constraint on 
the maximum possible value of a and fitting the _R-band 
light curve in two time portions 300 s < ti < 1000 s and 
t 2 > 800 s, we derive a>i < 0.88 ± 0.2 and a 2 = 0.31 ± 0.12. 
The second decay ai is found to be shallow whether or not 
the late time VLT i?-band measurement at 12.68 days is in- 
cluded and, by assuming that a break has not occurred in 
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Table 5. X-ray spectral indices for GRB 060108. Errors are at 90% confidence lev el. W hen possib le, we repo rt the 90% upper limit on 
the excess in column density, over the Galactic value of Ng al = 1.7 X 10 20 cm- 2 iPickev fc LockmanlH99oT) . The energy index of the 
first and second segment have been computed twice: with (3<) and without (ft) the inclusion of the X-ray flare. In the second case Cash 
statistics has been used. 



Segment 


Interval 


ft 


Nh excess 

_9 

cm 




1 

Flare 

2 

3 


t- ^ ^break,l 

^brcak,l ^ t ^ ^break,2 
t > tbreak,2 


1.31 ± 0.44 

0.68 ± 0.23 
0.71 ± 0.15 


< 1.4 X 10 21 

< 1.7 x 10 21 


1.47±°;« 
1.7 ± 0.4 
0.89 ± 0.3 



100 



Bessell B (FTN) 
Bessell V (FTN) 
Bessell R (FTN/VLT) 
SDSS i' (FTN) 
K (UKIRT/VLT) 
XRT (0.3-10 keV) 



100 1000 10000 100000 1e+06 

Time since GRB (s) 

Figure 4. BVRi' multicolour light curves of GRB 060108 ob- 
tained by the FTN pl us i\-ba nd detections from UKIRT (value 
taken from Levan ct al J2006rl and the VLT. Flux densities have 
been corrected for Galactic extinction. The XRT light curve is 
also shown for comparison. 



the optical light curve, this may represent a faint detection 
of the afterglow. On the other hand, as suggested by the 
X-ray light curve, it is possible that a break has already oc- 
curred around tbrcak,2 ~ 11-5 ks, in which case the predicted 
magnitude of the optical counterpart would be significantly 
fainter than R ~ 23.5 mag at 12.68 days. As discussed in 
§ 0] if taken at face value the apparent rebrightening in the 
i'-band light curve (Figure |1J may be interpreted as due 
to a further emission mechanism (possibly a second reverse 
shock) . In such a case, the decay at the end of the first hour 
may be significantly steeper than oli ~0.31, similar instead 
to ol\. Moreover, in the late VLT image the object is not 
pointlike but clearly extended, making likely that the emis- 
sion detected at such late time is dominated by the host 
galaxy (at least ~ 50% of the total flux). 



3.4-1 Photometric Redshift Derivation 

In order to estimate the photometric redshift of 
GRB 060108, we adopted a \ 2 minimization of the observed 
optical to infrared Spectral Energy Dist ribution (SED) 
based on the technique presented bv lFontana et alJ J2000). 
This is a widely used and well tested technique for redshift 
determination of galaxies (in which case it takes into ac- 
count for st ar formation history of each galaxy type and 
other effects: iFernandez-Soto et alJll999l . lCsabai et al . 2000, 
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Figure 5. Best fitting power law curves for the afterglow as ob- 
served with the FTN, consistent with the gradient between the 2 
NIR datapoints (not shown). Upper limits (including those regis- 
tered by the UVOT) are indicated with arrows. 



lRowan-Robinsonl200ilMelandri et all2006l) which has been 
recently readapted to fit simple power-law SEDs, as those 
typical of GRB afterglows. 

The photometric redshift was determined by combining 
contemporaneous optical and NIR flux measurements de- 
rived from the early-time FTN and UKIRT imaging at 45 
minutes after the burst. The BRi'K values were extrapo- 
lated to 3600 s using different values for the temporal decay 
index ranging from a =0.3 to a=0.6 (see F igure Et. Follow- 
ing the method described in lMelandri et alJ i2006tl . the mag- 
nitudes were converted to the AB system before inputting 
them into the photometric redshift determination algorithm. 

Adopting a \ 2 minimization of the observed spectral 
energy distribution (which is assumed to follow a power law 
F„ ex v" 13 ) at 3.6 ks after the trigger, we derive a maximum 
redshift of z < 3.2 at 90% confidence level. 

Moreover, by using the temporal decay value of a — 0.4 
obtained by the simple power law fit discussed in § 13.41 the 
minimization returns (3o ~ 1.42. The results obtained by 
using this value for a are shown in Figure EJ Here the top 
panel shows the derived AB magnitudes, corresponding to 
BRi'K values, and the best-fitting curve; the bottom panel 
shows the reduced \ 2 distribution as a function of photomet- 
ric redshift. The \ 2 distribution does not show a minimum: 
it is very flat from z ~ to z ~ 2 and then start increasing 
if z > 2.4. Therefore, it does only allow setting an upper 
limit determination for z. 
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Figure 6. Top panel: the optical:NIR spectral energy distribu- 
tion for contemporaneous magnitudes in B, R, i' and K bands. 
Observed values were extrapolated to 3600 s after the burst by 
using a decay slope of a = 0.4, coincident with that of the optical 
temporal decay. Bottom panel: the reduced \ 2 distribution as a 
function of photometric redshift. See text for details. 



The numerical technique used here does not account 
for intrinsic absorption. In the attempt to estimate this ef- 
fect, we imposed a power-law on the whole SED from the 
optical to X-ray wavelength range, and used a SMC extinc- 
tion profile. By taking z ~ 2, we find Ay 



Pox = 0.54 



+0.23 



0.26l»;?? and 



The shallower value of 0ox = 0.54 



+0.23 



is preferable: if [3o is as steep as ~ 1.1 we would be un- 
able to account for the SED from optical to X-ray accord- 
ing to the synchrotron spectrum of the afterglow. Moreover, 
Pox = 0.54lg'j2 i s mn y consistent with the X-ray slope 
registered at the same time, as shown in Table |K| 



4 DISCUSSION 

4.1 Afterglow Model 

The X-ray afterglow of GRB 060108 is quite complex 
and exhibits several of the interesting features which have 
been recently observed in others GRBs by Swift (see 



uift 

Zhan g et al .120061. iNousek et al]l2006l iBurrows et alJl2005bl 
Roming et al.l2006l for detailed discussions). In particular, in 
several cases the X-ray emission has been found to consist of 



an initial steep decay, followed by a shallow one about 100- 
1000 s after the trigger. A second change of the decay slope 
usually occurs about 10 ks later. Flares can be observed su- 
perimposed to this template behavior, and they are usually 
thought to be produced by internal shocks that should occur 
during the prompt emission phase. 

In the case of GRB 060108 we have registered three 
different spectral slopes. The X-ray lightcurve initially fol- 
lows the standard pattern of rapid to shallow transition, and 
then, after tb rC ak,2 ~ H-5 ks, the afterglow enters the "nor- 
mal afterglow" phase. Also, a small X-ray flare has been 
observed at t ~ 300 s. 

According to the standard interpretation, the initial fast 
decay is likely to be associated to the tail of the prompt emis- 
sion. In such a case, the relationship between the spectral 
and decay indices should be given by a = 2 + ("curva- 
ture relation", see iKumar fc Panaitescull2000l . Tzhang et alJ 
120061) . In the case of GRB 060108 it is ai = 2.78 ± 0.43 and 
0i — 1.47 ± 0.42, therefore the agreement is only marginal. 
The value of a required for the curvature relation to be sat- 
isfied is ~ 3.5, although the shallower initial slope observed 
in our case may possibly be caused by the presence of the 
flare at ~ 300 s. In fact, the a = 2 + relation gives the 
steepest decay index for a spherical model and is expected 
to be satisfied only when emission from on-axis fluid ele- 
ments is shut off suddenly, and off-axis emission dominates. 
The presence of minor internal shocks at later times would 
make the decay more gentle. On the other hand, Figure [3] 
shows that the backward extrapolation of the XRT data well 
matches the flux level of the prompt emission, as registered 
by BAT and extrapolated in the XRT band, strengthening 
the idea that the initial fast decay of the X-ray emission is 
related to the prompt emission. 

Just after the rapid decay phase, we observe a small 
flare that peaks at 300 ±10 s and, when modelled with a 
Gaussian, has a width oq = 32.3 ± 10.0 s. By adopting dif- 
ferent rebinning for the X-ray lightcurve, we have found that 
the x 2 always shows an improvement by adding a Gaussian 
to the fitting model, and that the confidence level for the 
flare is ~ 2 — 2.5<r. Such X- ray flares are common features 
of Sw ift XRT light curves dZhang et alJl200rl Eiang et alJ 
2006), and are usually associated with collisions between 
ultrarelativistic shells emitted by the "inner engine" after 
the main high-energy event. In this case, the decay slope 
and the spectral index of the flare should again satisfy the 
curvature relation, at least provided that the decay slope is 
computed by re-setting t he time-zero point to at the start 
of the shell emission (see IZhang et al.ll2006l and references 
therein). Although the flare statistic is quite low and does 
not allow a robust determination of the slope, we find that 
this relation is satisfied by taking to — 282 s. This gives 
otfi = 4.7 ± 1.1 (90% confidence level, while the spectral 
index is = 1.7 ± 0.4, see Table^. 

After tbrcak.i = 250 ± 30 s, the GRB 060108 afterglow 
enters a phase of plateau, which we can attribute to the 
emission from the forward shock produc ed when the fireball 
runs into the circu mburst medium (see IZhang et al]l2006l . 
lO'Brien et aljl2006t) . However, the observed temporal decay 
slope is shallow er than that exp ected from a standard af- 
terglow decline JSari et alJll998l) . requiring that additional 
energy is injected into the afterglow ("refreshed" afterglow). 
Three possible mechanisms have been proposed to explain 
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this late energy injection iZhang et al.ll2006l and references 
therein): i) a continuous activity from a longer lasting cen- 
tral engine, usually mimicked by a central engine luminosity 
law of the kind L oc t~ q , ii) a power-law distribution of the 
Lorentz factors in the ejecta, causing slower shells to catch 
up with the decelerated fireball or iii) the deceleration of 
highly magnetized ejecta, in which case the outflow is Poynt- 
ing flux dominated and transfers energy to the medium. In 
the latter case, however, a good parametrization of the en- 
ergy release is still lacking, theref ore we will not discuss this 
pos sibility in further d etails (see IZhang fc Kobavashil 1200*51 
and lZhang et al.ll2006l for a discussion). 

By using the closure relations listed by Zhang et al. 
(2006; see their Table 2), we find that the relation which 
is satisfied by the spectral and decay index during this 
phase (a ~ 0.3 and /3 ~ 0.7, see Tables [I] and is 
a = {q - 1) + (2 + q)P/2, giving q = 0.44. All other pos- 
sibilities are excluded. This implies that the fireball is prop- 
agating into a constant density medium (such as the inter- 
stellar medium, ISM) and that the X-ray frequency lies be- 
tween the synchrotron frequency and the cooling frequency 
(y m < v < i>c). During this epoch, the central engine does 
not turn off, but continues to inject energy in the outflow 
according to the law L oc t~ ' 44 . This result makes unlikely 
that the late energy injection is due to electromagnetic emis- 
sion from an highly magnetized and rapidly rotating pulsar 
(which in principle can be formed after the explosion of the 
supermassive star producing the GRB), as this would nor- 
mally require a flatter injection rate during the initial spind- 
down phas e (q ~ 0), followed by an asymptotic decay with 
q ~ 2 fsee lZhang fc MeszaroslfeoOll) . Similarly, any mecha- 
nism causing a nearly constant input of energy (for example 
early phases of disk accretion, during which the accretion 
rate, M, can be nearly constant, etc ...) can be ruled out. 
However, continued infall of matter onto a newly formed 
black hole may occur with q < 1, since at late times M is 
expected to fall off like a power law. 

As mentioned before, an alternative and observation- 
ally indistinguishable scenario is one in which the central 
engine activity is brief, but the ejecta have a distribution of 
Lorentz factors. The fastest shells initiate the forward shock, 
decelerate, and are successively caught by the slowest shells 
producing internal shocks. The consequent addition of en- 
ergy in the blast-wave mitigates the deceleration and the 
afterglow decay rate. Assuming a Lorentz factor distribu- 
tion M(> 7) oc 7 s , where 7 is the shell Lorentz factor, the 
index s is related to q by s = —(10 — 7q)/(2 + q), which, in 
this case, gives s = —2.83. 

The X-ray light curve breaks steeply at 11.5 ks. We can 
exclude that this is due to the passage of the cooling fre- 
quency through the X-ray band at this time, since in such 
a case, when comparing parameters relative to the epochs 
before and after the time break, we would expect to ob- 
serve a relatively large change in the X-ray temporal slope 
of Aa = 0.25, accompanied by a similar change in the spec- 
tral slope of A/3 = 0.5. In the case of GRB 060108, the 
observed change in the decay slope is Aa = 0.76 ± 0.12, 
while the spectral change is only A/3 = 0.03 ± 0.27 at 90% 
confidence level. This means that this interpretation is not 
compatible with the data (less than ~ 3cr). 

More likely, the second break at £break,2 = 11.5 ±0.4 ks 
occurs when the additional energy injection ceases and the 



afterglow enters a "normal" decay phase. During this phase, 
temporal and spectral indices satisfy the closure relation a = 
3/3/2, which again confirms that the fireball is propagating in 
an ISM-like medium. This also requires the X-ray frequency 
to be below the cooling frequency during the entire duration 
of XRT observations. The corresponding index of the power- 
law energy distribution of radiating electrons is p = 2.39. 

The comparison between the optical and X-ray light 
curve is also illuminating. Interestingly, the i! and R band 
flux after the first 800 s have a decay slope of ~ 0.3 — 0.4, 
similar to what is observed at the same time in the X-ray 
band. This suggests a common origin for the X-ray and op- 
tical emissions at these epochs. 

On the other hand, detailed observations in the R and 
i' band may suggest that the optical decay slope may be 
steeper between 300 and 800 s. This may indicate that 
we are detecting a decay from an early optical "bump", 
i.e. that at such early times we are observing the fading 
emission from a faint reverse shock, which propagates back- 
ward into the ejecta w hile they are refres hed by the late 
central engine activity iZhang et al. | |20031. A possible re- 
brightening is observed in the i' band only. The signifi- 
cance is low and a monotonic decay is not ruled out. On 
the other hand, it is worth to point out that a similar fea- 
ture has been detected at larger signifi cance in the brighter 
GRB 060206 ( Monfardin i et al . 2006c), in which case the R 
and i' lightcurves could not be fitted with identical func- 
tions. If real, the rebrightening in the i' band may de due to 
a second reverse shock, possibly initiated by the collision of 
the 300 s of the X-ray flare into the external medium. 

The X-ray light curve breaks steeply at 11.5 ks, while 
no apparent corresponding break is observed in the i' and 
R bands. A time break produced by the passage of the cool- 
ing frequency through the X-ray spectrum would certainly 
account for this, since in such a case the break would be 
monochromatic. 

On the other hand, based on the spectral and temporal 
analysis, we find that this is not likely to be the case and the 
X-ray break is most probably due to the cessation of energy 
injection. In such a case, and if the optical and X-ray emis- 
sion were produced by the same mechanism, a break should 
be expected to appear in the optical emission, at ~ ibreak,2- 
Here we only notice that observations are sparse at this time; 
there are no data in the i' band measured after 5ks and 
the only _R-band measurement taken after 11.5 ks indicates 
that the light curve does not break. The fit of the optical 
lightcurve with a single power law suggests that an optical 
break may not be required even in correspondence of the first 
X-ray break. It is therefore possible that the GRB 060108 
behaviour is similar to that of other GRBs, which do not 
present a bre ak in the optical contem poraneous to that in 
the X-ray fsee lPanaitescu et al.ll2006T) . A satisfactory expla- 
nation for this behaviour is still lacking. Alternatively, if the 
i?-band observation (7?=23.5 mag taken at 12.67 days after 
the trigger) is attributed to emission from the host galaxy of 
the GRB, a break in the optical GRB light curve at ibreak,2 
would not have been observed. As mentioned before, a sig- 
nificant host galaxy contamination is expected, since the 
source appears to be extended in the late VLT image. 
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Figure 7. Comparison of GRB afterglow light curves with that 
of GRB 060108, illustrating the observed faintness of its optical 
emission. Re-adapted from Guidorzi et al. (2005); see the original 
paper for the references of all data. 

4.2 Observed Darkness 

Figure[7|shows the lightcurve of GRB 060108 in the broader 
context of other GRBs. One peculiar characteristic of this 
event is that, although the burst is detected in the BVRi' 
bands at later times, the UVOT does not observe any opti- 
cal/UV emission within the first ~ 300 s (i.e. during the tail 
of the prompt emission where internal shocks are expected 
to occur). 

In order to investigate the possible reason for this 
behaviour, we can compare the relative faintness of 
GRB 060108 with other bursts using the X-ray flux ver- 
sus gamma-ray fluence and the optical flux versus X-ray 
flux plots published by iRoming et alJ (120061) . In the case 
of GRB 060108, at 1 hr after the burst the X-ray flux 
is 3.1 x 10 _12 erg cm _2 s _1 and the gamma-ray fluence is 
5" = 4.9 x 10" 7 ergcm- 2 (in the 2-10 keV and 15-350 keV 
bands, respect i vely). As we can see from the Figure 3 of 
IRoming et al] <l2006l) . this burst is not unusual and sits 
solidly within 2a of the mean, although overall it is slightly 
fainter in both 7 and X-rays compared to other GRBs. We 
therefore expect a typical behaviour also in the optical emis- 
sion (perhaps with a slightly fainter optical afterglow than 
other bursts). If we compare the X-ray and optical fluxes at 
1 hour after the trigger (when the optical flux is ~ 8.0 ^Jy)) 
we find values within the range of the population with a 
slight X-ray to optical excess. Similarly, if we extrapolate the 
optical point at 11 hours, we find an optical flux of ~ 3.9 fj,Jy, 
while the X-ray flux is 4 x 10 -13 ergcm _2 s _1 . Although we 
caution that the behaviour of the light curve at 11 hr is not 
well constrained (see below), if taken face value this would 
imply that GRB 060108 still lies in the 2a limit, near the 
region of optically faint events. 

A further w ay to assess the "darkn ess" of a GRB has 
been proposed bv ljakobsson et all ll2004l) (see also lRol et al] 
2005), on the basis of the optical to X-ray spectral proper- 
ties. According to these authors, a burst can be classified 
as "truly" dark when Pox < 0.5. A ratio less than this 
threshold would be in disagreement with the spectral slopes 
expected from the synchrotron emission, requiring an intrin- 
sic "lack" of optical emission. Indeed, the ratio of X-ray and 
R band measurements of GRB 060108 provide a value of 



Pox = 0.45, as measured at ~ 100 s afte r the burst trigger. 
According to the Jakob sson et al] (I2004T) classification, this 
burst is therefore "dark" (see Fig. |5J . 

In principle, a possible explan ation for the "observed 
darkness" is a large redshift (z > 5.lLamb fc ReichartfeOOCl 
iKawai et al ] l2006t iGendre et al]l200rj) . In such a case the 
UV light, which is heavily attenuated by intergalactic hy- 
drogen, is redshifted into the optical band. The burst would 
then appear very red. However, the high redshift scenario 
is ruled out in this case as we observe the afterglow in 
the B band, and we have been able to determine an up- 
per limit of the redshift of z < 3.2 by using the pho- 
tometric redshift technique. High magnetic fields in the 
ejecta may also restrict the amount of energy transferred to 
electrons by reverse shock s, inhibiting the optical emission 
JZhang fc Kobavashil20(55l) . 

Alternatively, the burst could be obscured by dust. Dust 
would cause extinction of the afterglow and the emission 
appears reddened as the shorter wavelengths are absorbed 
more. It is interesting to compare the photoelectric absorp- 
tion measured from the X-ray spectrum with the reddening 
inferred from the optical. As discussed in H3.3I by using all 
data after the flare through to the end of the observations, 
we have been able to measure a significant column density 
excess (99% confidence level). Assuming z = 2, this yields 
N H = 8.6l 7 ; 2 x 10 21 cm" 2 . 

Light detected in the R band, centered on 700 nm in 
the Earth restframe, is emitted at 233 nm in the burst lo- 
cal frame (by assuming again z ~ 2). Assuming a SMC 
extinction curve, and Nn/A y — 1.6 x 10 22 cm~ 2 /mag 
( Wcingart ner fc Draind |2J}0Qj) , the extinction determined 
from the optical-IR SED in Section 3 (Ay = 0.26±g 5? ) cor- 
responds to an X-ray derived photoelectric absorbing col- 
umn density of N H = x 10 21 cm 2 consistent with 
the X-ray measurement. Note however, that for the SMC 
metallicity, considerably lower X-ray absorption would be 
expected for this Ay - This amount of extinction is sufficient 
to bring the optical emission up to a value compatible with 
regular Swift bursts. Dust extinction would also alter the 
spectral index of the optical decay. Our estimated value, if 
accounted for, would bring the Jakobsson flux ratio above 
0.5, thus removing this burst from the optically dark regime. 

The possibility that the "observed optical darkness" is 
caused by intervening dust along the line of sight is also 
favoured by results from the optical - IR SED modeling. 
As discussed in § 13.4.11 without accounting for extinction 
the technique gives a steep Po = 1.42, that would be hard 
to explain in terms of a synchrotron spectrum of the af- 
terglow. When extinction is included, the value of Pox is 
0.541q 12) which compares better with the observed SED. It 
can be noticed that the value of Pox is already low, close to 
synchrotron emission limit of 0.5, even after correcting for 
extinction. The hardness of the spectrum is another factor 
that contributes to make this burst optically faint. 

It would then appear from the arguments above that a 
likely explanation for the observed darkness is a combination 
of intrinsic faintness, hard optical to X-ray spectrum and 
dust extinction. 
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Figure 8. Ratio of optical to X-ray brightness as a function 
of time of GRB 060108, compared with other GRBs, dark and 
luminous. Re-adapted from Jakobsson et al. (2004). 



5 CONCLUSIONS 

We presented and discussed the gamma-ray, X-ray and 
UV /optical/infrared observations of the Swift GRB 060108, 
performed with the instruments on-board the spacecraft as 
well as the with the ground-based Faulkes Telescope North, 
United Kingdom Infrared Telescope and Very Large Tele- 
scope. GRB 060108 has a moderately faint X-ray afterglow, 
but very low optical emission, making it one of the "dark- 
est" GRBs ever observed. We suggest that faintness of the 
optical emission may be due to an intrinsic weakness of the 
burst and a hard optical to X-ray spectrum, accompanied 
by some degree of extinction which occurred in the GRB 
surroundings. 

The X-ray light curve shows the typical template 
discovered by Swift, characterized by a rapid decay in the 
first 300 s (commonly interpreted as the tail of the prompt 
emission on the basis of the combined temporal/spectral 
properties) followed by a fiat decay slope (typical of "re- 
freshed" afterglows). An X-ray flare is observed at ~ 300 s, 
although the statistic is low. The flat decay phase lasts 
for about ~10 ks, after that the light curve breaks into 
a steeper segment with a power-law decay index of 1.0, 
typical of a standard afterglow phase. 

A comparison with the optical light curve is quite inter- 
esting. While the optical flux after 800 seconds has a slope 
of a = 0.4, similar to that of the X-ray in the same inter- 
val, the decay index before this time is likely to be steeper. 
This behavior may be explained if we interpret the fast de- 
cay as the reverse shock emission, while the later, flatter 
emission, has the same origin of the X-ray one. A possi- 
ble rebrightening is visible only in the infrared, although 
a monotonic decay is not ruled out. A si milar feature has 
been seen in the case of GRB 060206 (iMonfardini et alJ 
l2006d) . If real, it might indicate a second reverse shock 
emission, initiated by the 300 s flare. Alternative models 
able to explain an optical rebrightening include an increase 
in the densit y of the medium w here the forward shock 
is produce d jLazzatLeLal] |2002|) , or energy injection by 
late shells llBior nsson 1 Gudmundsson fc J6hannessonll2004l 
I J6hanne5scmet^rtTl27)0a) . This last scenario is less likely, be- 
cause of our finding that the late energy injection is taking 
place at a steady rate. The coincidence of the rebrighten- 



ing with the X-ray flare, however, favours the hypothesis of 
reverse shock. 

Another intriguing feature of the optical light curve is 
the absence of a break in correspondence to the X-ray one at 
~10ks, which may be either intrinsic or due to a significant 
contribution to the optical flux at late times by the host 
galaxy. However, the poor sampling of the late optical light 
curve does not allow us to better constrain this behaviour. 

The analysis of the optical spectrum, obtained from 
data gathered 45 minutes after the burst, has allowed us 
to determine an upper limit of z < 3.2 at 90% confidence 
level, by using a \ 2 minimization of the observed spectral 
energy distribution. 

The optical afterglow is below the detection limit of the 
UVOT within 100 s of the burst. This event has shown how 
observations taken promptly and deeply enough may reveal 
the interesting behaviour of the early optical emission. In the 
Swift era, further similar observations are a reality due to the 
prompt response of both the spacecraft and ground robotic 
telescopes, as well as due to the possibility to perform deep 
observations with large telescopes at early times after the 
trigger. 



6 ACKNOWLEDGEMENTS 

CGM acknowledges financial support from the Royal Soci- 
ety. SZ thanks PPARC for support through an Advanced 
Fellowship. KP thanks PPARC for support. AM acknowl- 
edges financial support from the Provincia Autonoma di 
Trento. The Faulkes Telescope is operated with support from 
the Dill Faulkes Educational Trust. The UKIRT is operated 
by the Joint Astronomy Centre on behalf of the UK PPARC. 
VLT observations were carried out under programs 076. A- 
0392 and 076.D-0667. We acknowledge excellent support 
from the observing staff in Paranal. We thank the referee, 
Prof Ralph Wijers, for useful comments and suggestions. 



REFERENCES 

Akerlof, C. et al., 1999, Nature, 398, 400 

Barthelmy S.D. et al., 2005, ApJ, 559, 710 

Bessel, M.S., 1979, PASP, 91, 589 

Barthelmy, S.D. et al. 2005, Nature, 

Berger, E. et al., 2002, ApJ, 581, 981 

Berger, E. et al., 2005, ApJ, 634, 501 

Bertin, E., Arnouts, S., 1996, A&AS, 117, 393B 

Bjornsson, G., Gudmundsson, E.H., & Johannesson, G., 

2004, ApJ, 615, L77 
Bloom, J.S. et al., 2006, GCN Circ. 4475 
Blustin, A.J. et al, 2006, GCN Circ. 4450 
Bromm, V. & Loeb, A., 2002, ApJ, 575, 111 
Burrows D.N. et al., 2005, Space Science Review, 120, 165 
Burrows D.N., Romano P., Godet O. et al., 2005b, 

in press in the proceedings of "X-ray Universe 2005" 

(astro-ph/0511039) 
Butler, N. & Bloom, J.S., 2006, GCN Circ. 4500 
Cardelli et al. 1989 

Cool R.J. et al., 2006, GCN Circ. 4447 
Cool R.J. et al., 2006, AJ, 118, 733 
Costa,E. et al., 1997, Nature, 387, 783 



12 Oates et al. 



Covino, S. et al., 2006, A&A, in press 

Csabai I., et al., 2000, AJ 119, 69 

Dai Z.G., Lu T., 1998a, A&A, 333, L87 

D'Avanzo, P., et al., 2006, GCN Circ. 4501 

De Pasquale, M. et al., 2003, ApJ, 592, 1018 

De Pasquale, M et al. 2006, MNRAS, 365, 1031 

Dickey, J.M., Lockman, F.J., 1990, A&A, 28, 215 

Fernandez-Soto A., et al., 1999, ApJ 513, 34 

Fioc M. & Rocca-Volmerange B., 1997, A&A 326, 950 

Fontana A., et al., 2000, AJ, 120, 2206 

Fox, D.B. et al., 2005, Nature, 437, 845 

Fruchter, A.S, 1999, ApJ, 512, LI 

Fukugita, M., Ichikawa, T., Gunn, J.E., Doi, M., Shi- 
masaku, K., Schneider, D.P, 1996, AJ, 111, 1748 

Fynbo, J.U. et al, 2001, A&A, 369, 373 

Gehrels, N. et al., 2004, ApJ, 611, 1005 

Gehrels, N. et al., 2005, Nature, 437, 851 

Gendre, B. et al., 2006, A&A, submitted, 
(astro-ph/0603431) 

Groot, P.J. et al., 1998, ApJ, 502, L123 

Guidorzi, C. et al., 2005, ApJ, 630, L121 

Guidorzi, C. et al., 2006a, PASP, 118, 288 

Guidorzi, C. et al., 2006b, GCN Circ, 4447 

Jakobsson, P., Hjorth J., Fynbo J.P.U. et al., 2004, ApJ, 
617, L21 

Johannesson, G. et al., 2006, ApJ accepted 

(astro-ph/0605299) 
Kalberla P.M.W., Burton, W.B., Hartmann, D., et al. 2005, 

A&A, 440, 775 
Kawai, N. et al., 2006, Nature, 440, 184 
Kumar, P., & Panaitescu, A., 2000, ApJ, 541, L51 
Lamb, D.Q. & Reichart, D.E., 2000, ApJ, 536, 1 
Lamb, D.Q. et al., 2004, New Astronomy Reviews, 48, 423L 
Lazzati, D., et al. 2002, A&A, 396, L5 
Lazzati, D., Covino, S. & Ghisellini, G., 2002, MNRAS, 

330, 583 

Le Borgne D. & Rocca-Volmerange B., 2002, A&A 386, 446 

Levan, A. et al. 2006a, GCN Circ. 4503 

Levan, A. et al. 2006b, GCN Circ. 4457 

Liang, E.W., Zhang, B., O'Brien, P.T. et al. 2006, ApJ in 

press (astro-ph/0602142). 
Melandri, A., e t al., 2006, A&A in press, 

(astro-ph/0510604) 
Monfardini, A. et al., 2006a, GCN Circ. 4502 
Monfardini, A., et al. 2006b, GCN Circ. 4630 
Monfardini, A., Kobayashi, S., Guidorzi, C, et al. 2006, 

ApJ in press, ( astro-ph/0603181 ) 

Nousek J. et al., 2006, ApJ, 642, 389 

O'Brien P.T. et al., 2006, ApJ, in press, ( astro-ph/0601125 ) 

Oates S.R. et al., 2006, GCN Circ. 4443 

Page K.L. et al., 2006, GCN Circ. 4453 

Panaitescu, A., Meszaros, P., Burrows, D., Nousek, J., 

Gehrels, N., O'Brien, P., Willingale, R., 2006, MNRAS, 

369, 2059 
Pedersen et al., ApJ, 636, 381 
Piranomonte S. et al. 2006, GCN Circ. 4484 
Poole T.S., et al., ApJ, 636, 381 
Rees M.J., Meszaros P., 1998, ApJ, 496, LI 
Reichart D.E., Price. P.A., 2002, ApJ, 565, 174 
Rol, E., Wijers, R.A.M.J., Kouveliotou, C. et al., 2005, 

ApJ, 624, 868 
Roming P. et al., 2006, |astro-p h/0509273 



Rowan-Robinson M., 2003, MNRAS 345, 819 
Sari R., Piran T., Narayan R., 1998, ApJL, 497, 17. 
Schlegel, M. F., et al., 1998 
Skrutskie, M. F., et al., 2006 AJ 131, 1163S 
Smith, J. A. et al., 2002, A J 123, 2121 
Steele, I.A., 2004, Astronomische Nachrichten, 325, 519 
Tagliaferri, G. et al., 2005, A&A, 443, LI 
van Paradijs, J. et al., 1997, Nature, 386, 686 
Villasenor, J.S. et al., 2005, Nature, 437, 855 
Weingartner J.C. & Draine B.T., 2001, ApJ, 548, 296 
Zhang B., Meszaros P., 2001, ApJ, 552, L35 
Zhang B., Kobayashi S. & Meszaros P., 2003, ApJ, 592, 
950 

Zhang B., Meszaros P., 2004, IJMPA, 19, 15, 2385 

Zhang B., Kobayashi S., 2005, ApJ, 628, 315 

Zhang B., Fan Y.Z., Dyks J. et al., 2006, ApJ, 642, 354 



